Creation, manipulation and detection of Majorana fermions with cold atoms in optical 

lattice 
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We propose an experimental scheme to simulate the transverse field Ising model with cold atoms 
trapped in one-dimensional optical lattice. Majorana fermions are created at the ends of the optical 
lattice segment in topological phase. By controlling the addressing lasers, one can move, fuse and 
braid them. We also show that the non-Abelian braiding statistics of Majorana fermions can be 
demonstrated unambiguously through the construction of two braiding operations and distinguishing 
the resulting two output orthogonal collective spin states. A nice feature of the scheme is that the 
strong fluorescence provided by the collective spin state can be readily detected in experiment. 
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Majorana fermions are particles that are their own 
antiparticlcs unlike Dirac fermions where electrons and 
positrons are distinct [l| . When exchanged among them- 
selves, Majorana fermions obey non-Abelian statistics. 
Such fundamentally interesting particles have recently 
attracted much attention due to their potential appli- 
cation in topological quantum computation Q. Majo- 
rana modes are originally perceived as zero-energy states 
bound to the vortices in two-dimensional (2D) spinlcss 
p x ± ipy-WAve superconductor (SPSC) [3|l or the two ends 
in one-dimensional (ID) SPSC chain [J]. It has been 
recently proposed that a semiconductor thin film with 
Rashba spin-orbit coupling, together with proximity- 
induced superconductivity and Zceman splitting, resem- 
bling the SPSC model, can be used to create Majorana 
fermions [f| . One dimensional version of this system has 
been shown to host Majorana fermions at the two ends 
of a semiconducting wire Q. By using a T-junction wire 
network, the Majorana fermions in the ID semiconduct- 
ing wires can be braided by tuning the local gates Q- 
Despite the relative ease with which Majorana fermions 
can be stabilized in ID wires making this ID set-up more 
promising as a topological quantum information process- 
ing platform, experimental challenges need to be over- 
come: strong Zeeman fields could destroy the supercon- 
ductor and the local control of electron density by the 
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gates is non-trival due to strong screening by the super- 
conductor. 

Cold atoms nowadays are widely recognized as pow- 
erful experimental tools for mimicking a wide range 
of systems originally stemming from condensed-matter 
physics. This system can provide a highly controllable 
and tunable environment. Many significant experimental 
advances have been made in this forefront field. In par- 
ticular, recent experiments have realized synthetic mag- 
netic fields and spin-orbit coupling for ultracold atoms 
[H E|. Several protocols have been put forward to use 
these technologies for generating and pro bing Majorana 
fermions in 2D and ID fermionic gases p^J - [l5| . However, 
the braiding Majorana fermions as well as the detection 
of their non-Abelian statistics remains an outstanding 
problem. 

In this work, we propose an experimental scheme to 
create and braid Majorana fermions as well as to de- 
tect their non-Abelian statistics with cold atoms in opti- 
cal lattice. Note that the transverse field Ising model 
(TFIM) after Jordan- Winger (JW) transformation is 
equivalent to the ID SPSC model [II, We show that 
various TFIMs can be simulated with two-component ul- 
tracold bosonic atoms trapped in a ID optical lattice. 
This atom-lattice system is highly tunable and very ro- 
bust. The spin exchange coupling and the transverse 
field can be tuned by changing the laser intensity. Un- 
der this tuning, the optical lattice can be driven into 
topological phase or non-topological phase. Majorana 
fermions are bound to the ends of the optical lattice seg- 
ment in topological phase. Such tuning can also allow 
Majorana fermions to be moved, recreated and fused. 
By further employing a cross lattice, Majorana fermions 
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can be braided in the ID optical lattice. In addition, the 
orientation of our simulated Ising coupling can be conve- 
niently rotated by varying the phase of the laser. Finally, 
we show that such rotation can also mimic the braiding 
of Majorana fermions. 

Unambiguous detection of non-Abelian statistics of 
Majorana fermions through easily measurable collective 
spin states is another distinct advantage of our proposal. 
We use the braid group elements of Majorana fermions 
to construct two different orders of braiding operations. 
When applied to the same initial state, the correspond- 
ing output states are orthogonal. For ID superconductor 
model, the state of the non-local Majorana fermions has 
neutral charge and it is hard to be detected. However, 
the Majorana fermion states in our spin chain model can 
be mapped back to spin basis so that the states of two 
output Majorana fermions are two orthogonal collective 
spin states. Compared with the detection of the fluores- 
cence from single atom [it! [Tlj . which would be hard, 
these collective states consist of hundreds of atoms and 
they can provide strong fluorescence, which allows for ex- 
perimental detection with high efficiency. In this way, our 
proposed scheme provides an easy way to detect the fun- 
damental non-Abelian statistics of Majorana fermions. 

Let us consider an ensemble of ultracold Rb atoms 
trapped in a ID optical lattice at 1064 nm. With stan- 
dard optical pumping methods, we suppose that the 
atoms can only be in the states \5Si/ 2 ,F = 2, mp = 1) 
or \5Si/ 2 ,F = l,mp = 0) denoted in the following text 
with an effective spin index a =4,, j. Taking into account 
the ground state hyperfine splitting Auhf = 6.8 GHz, the 
AC polarizabilities for the spin up and spin down have a 
small difference of about 3%. As shown in Fig. 1, it can 
be used to generate a state-independent optical lattice by 
applying a standing wave laser beam L\. By considering 
a trapping depth of 20E r , where E r is the lattice re- 
coil energy, the ID potential is almost state-independent 
within a mismatch on the hopping rate lower than 15%. 
In order to simulate the Ising model in this ID optical 
lattice, we need to apply a second state-dependent op- 
tical lattice with a potential V a . As we will see, this 
state-dependent potential is very important for realizing 
the Ising model. 

At sufficiently low temperatures and strong lattice po- 
tential, the atoms are confined to the ground state of 
the optical lattice, and the system can be described by a 
Bose-Hubbard Hamiltonian 
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where ai a is the bosonic annihilation operator for atom 
states of spin a at the lattice site i, rii a = af a ai a , 
: (...): denotes normal order of the product of creation- 
annihilation operators, t a is the spin-dependent hopping 
rate and U aa i is the on-site energy between species a and 
a' . It is known that t a depends sensitively upon the lat- 
tice potential while U aa i exhibits weak dependence [l8j |. 
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FIG. 1: (Color online) Energy diagram and laser beam con- 
figurations for the a z (a) and (b) Ising models. L\ and 
L2 are the lasers leading to the ID state-independent optical 
lattice and spin-dependent hopping rate, (c) 2D lattice for 
the braiding of Majorana fermion. The filled blue (open red) 
circles symbolized the atoms in the \o) (|x)) state. An RF- 
field is supposed to be used to flip the state in \x) to the spin 
state \<j). 



Moreover, for the state under consideration the scatter- 
ing lengths are Off ~ 044. — . Thus, in what follows, 
U atJ ' = U is assumed. By varying the intensity of trap- 
ping laser to change the lattice potential, one can control 
the ratio U /t a to realize the Mott insulator with unitary 
filling [l9| . we will now consider this case in our model. 

By utilizing an additional statc-dcpcndcnt potential, 
the above Hubbard Hamiltonian is reduced to the a z - 
and (T^-Ising models. For <r z -Ising model, as sketched 
in Fig. 1(a), we apply a weak standing- wave laser beam 
L2, detuned between the D2 (780 nm) and Dl (794 nm) 
lines. In order to make complete overlap with the first 
state-independent optical lattice, the standing-wave laser 
beam should have an proper angle of incidence to the first 
one. If one tunes the intensity of the laser beam so that 
Vf S> Vi (or vice- versa), the hopping rate becomes 
negligible while t± remains finite, after compensating ef- 
fective B-fields, the effective Hamiltonian of Eq. (1) could 
be reduced to the crMsing model J\ J2j=i X a j a j+i w ^ n 
Ji = —t\/2U [l8j]. This method requires an interferomet- 
ric control of the spatial overlapping of the two standing 
waves at the atomic cloud position, which may be exper- 
imentally challenging. 

An alternative and simpler method, for cr^Tsing 
model, consists of coupling the effective spin states via a 
two-photon Raman process. One photon is issued from 
the ID lattice beams with a Rabi frequency Vl\ cos 2 kx, 
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where k is the wave number of the laser. The sec- 
ond photon comes from a running wave, L 2 , propagat- 
ing along x > with a complex Rabi frequency f/be 1 ^ 
(see Fig. 1(b)). At the equilibrium position of the 
atoms, the complex Rabi frequency of the Raman cou- 
pling is £Ir = ili^e 1 ^ I A, where A is the one photon 
detuning. One notes that ip can be controlled by stan- 
dard phase lock techniques and adjusted over a tempo- 
ral delay of (47rAi//i/) -1 ~ 12 ps. At resonance, the 
Raman coupling generates the state-dependent dressed 
potential V 2 = V+ |+) (+| + V_ |-) (-|, where |±) = 
(It) ± \i))/V2, V+ = -V_ = h\n R \/2. The to- 
tal state-dependent potential barriers now become V a = 
Vq + V s (s = ±). By improving the intensity of laser beam 
L-2 to make V+ S> VI, the hopping rate t+ for the atom 
in the dressed state |+) becomes negligible. This can be 
achieved with moderate laser power. In this case, the 
effective Hamiltonian of Eq. (1) can be rewritten as the 

cr^-Ising model J 2 Ylj=i (T j a 'j+i> where J 2 = —t 2 _/2U 
and cr v = cos (pcr x + sin <pa y . 

To realize the TFIM, we still need to generate the 
transverse field. For transverse field erasing model, an 
addressing laser beam with Rabi frequency Q x is applied 
to the lattice sites, which can generate the transverse 
field term h x cr x . Combined with the Ising Hamiltonian, 
we realize the transverse field cr z -Ising model 

JV-l N 

a : • ( 2 ) 

3=1 3=1 

where the transverse field h x = Ml x /2. For transverse 
field c^-Ising model, the transverse magnetic field h z <j* 
can be generated by applying an external magnetic field 
or by detuning the 2 photons Raman. We get a state- 
dependent AC shift V z <jj, which is equivalent to the ef- 
fective transverse magnetic field. Thus we arrive at the 
transverse field a^-Ising model 

JV-l N 

ff 2 = J 2 £«+l +/l *5>3*> (3) 
3=1 3=1 

where the transverse magnetic field h z = V z . In partic- 
ular, when the laser phase <p is tuned to or n/2, the 
above model is with respect to the transverse field a x - 
or (T^-Ising model. As one can see, our simulated TFIM 
is highly tunable based on controlling the lasers, includ- 
ing tuning the orientations of Ising interaction, the spin 
exchange coupling and the transverse magnetic field. 

Kitaev has shown that the transverse field cr^-Ising 
model after a JW transformation is equivalent to the 
ID SPSC model with superconducting phase = 
This result also holds true for transverse field cr v - and 
(T z -Ising models. Here we further show that our sim- 
ulated transverse field cr^-Ising model can be reduced 
to ID SPSC model with tunable superconducting phase. 



For our purpose, through employing the string-like an- 
nihilation and creation operator aj — aj Yii=i °f an d 

a+ = cr+ Ilti °f i tne transverse field cr v -Ising model 
can be rewritten as 

JV-l 

H 2 = J 2 ^(e l %-e^a+)(e^a J+ i+e-^a+ +1 ) 
3=1 

JV 

+fc,£(2at 0j --l), (4) 
3=1 

where the Dirac fermions satisfy the anticommutation 
relationship {a^a" 1 "} = Sij. By associating the spin ex- 
change coupling J 2 to the hopping amplitude and the 
superconducting gap, the laser phase tp to the supercon- 
ducting phase and the magnetic field h z to the chemical 
potential, the transverse field er^-Ising model is mapped 
into the ID SPSC model with tunable superconducting 
phase. Based on the N Dirac fermions annihilation and 
creation operators, 2N Majorana fermions operators are 
defined as 

jA,j = e l<p aj +e~ l(p af, j B ,j = i{e^ l<p af - e ilp aj), (5) 

where 7^ = 1 and {j a ,j, 7^} = 25 a/3 5 jk , a, /3 = A, B. 
In terms of these operators, the Hamiltonian in Eq. (4) 
becomes 

JV-l JV 
H 2 = 1J2 2J lB,jlA,j + l + ih z ^jA,jjB,j- (6) 
3 = 1 3 = 1 

When J 2 7^ and h z = 0, the system is in the topo- 
logical phase with two unpaired Majorana end modes 
7A,i and Jb,n, while for J 2 = and h z ^ 0, there is 
no Majorana end modes and the lattice is in the non- 
topological phase. In the previous case, the two Majo- 
rana fermions can be combined into an ordinary non-local 
Dirac fermion c = (7,4.1 + ^1b,n)/2, which yields two de- 
generate ground states |0) and |1) = c + |0) 0]. Because 
the fermion parity P = — ifl^Lj. lA.jlB.j has a Z 2 sym- 
metry, the ground states |0) and |1) have even and odd 
parity, i.e. P |0) (|1» = |0) (- |1)). Due to the nonlocal- 
ity of Majorana fermions that protects the ground states 
from decoherence, one can use them to encode a topolog- 
ical qubit for topological quantum memory. 

The spin exchange coupling J 2 and the transverse 
field h z can be tuned by controlling the addressing laser 
beams. As displayed in Fig. 2(a-b), such tuning allows 
Majorana fermions to be created, moved and fused in 
the ID optical lattice. In Fig. 2(a), we assume that the 
lattice is a simulated transverse field cr^-Ising model and 
we can address its two edge sites with two laser beams, 
one for improving the local state-independent potential 
to make J 2 = and one for generating the transverse 
field h z . The two edge lattice segments are driven into 
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FIG. 2: (Color online) Schematic setup for manipulating and 
braiding Majorana fermions. By controlling the addressing 
lasers, one pair of Majorana fermions (71, 72) is created and 
moved in (a) and (b). With the same control, two optical 
lattice segments L and R are driven into topological phase 
and two pairs of Majorana fermions (71, 72) and (73, 74) are 
generated in (c). To braid these Majorana fermions, a cross 
lattice is introduced in (d), where the orange segment de- 
notes the topological phase while the blue one denotes non- 
topological. 



non-topological phase but the center segment remains in 
the topological phase, which leads to one pair of Ma- 
jorana fermions 71 and 72 created at the ends of the 
topological segments. In Fig. 2(b), by sequentially driv- 
ing the leftmost lattice segment non-topological, the Ma- 
jorana fcrmion 71 is moved rightward and finally fused 
when it meets the Majorana fermion 72. The experi- 
mental requirement of single site addressing is now pos- 
sible with high resolution microscope [2(| [H( . As we will 
see later, this manipulation constitutes basic operations 
needed to exchange Majorana fermions and realize their 
non-Abclian braiding statistics. 

Another feature of our model is that the supercon- 
ducting phase for JW fermions is tunable. Interest- 
ingly, we find that this tuning can be used to simu- 
late the braiding of the Majorana fermions at the ends 
of the same topological segment. As described in Fig. 
2(a), two Majorana fermions ^a,i an d Jb,n arc lo- 
cated at the ends of the center topological lattice seg- 
ment. According to Eq. (5), they can be written as 
7A,i = e^c + e-^c+, -/ b ,n = i{e~ iip c + - e^c). By 
adiabatically tuning the laser phase from tp to tp ± 7r/2, 
we get the following transformation Ja,i ~ * T7b,jv, 
1b, n ±7i4,ij which just implements the braiding op- 
eration t = exp(=Ff Jb,n1a,i) = exp(q=f a z ) 0, where 
z is a Pauli matrix in the fermionic occupation basis 
{|0), c + |0)}. Moreover, adiabatic change of the laser 
phase tp by tt will result in a a z quantum gate operation. 
These operations are expected to be used for manipu- 
lating the quantum information stored in the topological 
quantum memory. 

We now consider the detection of the non-Abelian 
braiding statistics of Majorana fermions. It is known 
that, if one exchanges the Majorana fermion 7^ with 



its nearest neighbor one 7^+1, one can realize the braid- 
ing operation Tj = exp( : |7 i+1 7 i ) [22J. The non-Abelian 
braiding statistics arises from the fact that TiTi+\ 7^ 
T,-|_iTj. So in order to demonstrate this non-Abelian 
braiding statistics, we at least need four Majorana 
fermions. As shown in Fig. 2(c), by addressing some 
particular optical lattice sites, the two optical lattice seg- 
ments L and R are driven into topological phases and 
two pairs of Majorana fermions (71,72) and (73,74) are 
created at their ends. Based on these Majorana end 
modes, we define two Dirac fermions = (71 + 1^2)/^ 
and cr = (73 + 174 )/2. In Fig^ 2(d), as the way in 
the superconducting wire model [7[ , we introduce a cross 
lattice and use the basic operations in Fig. 2(a-b) to 
braid the Majorana fermions at the ends of the same 
or different topological segments, which correspond to 
the braiding operation s tu x) and t%. Following the ex- 
perimental method in [231 ] . as shown in Fig. 1(c), the 
braiding operations involving two orthogonal optical lat- 
tice segments can be realized by using a 2D Mott in- 
sulator with unitary filling. Initially the atoms in the 
2D Mott state are in a spin state \x) 7^ Using 
single spin addressing, the two orthogonal lattice seg- 
ments are created in the \a). If V x ,o Va,a'i the dy- 
namic is frozen between \x) and \a) neighboring sites. 
Here we assume the lattice has simulated the trans- 
verse field cr x -Ising model. In the fermionic occupation 
basis {|00) ifl) c\ |00) Li? , c+ |00) ii? , c+c+ |00) iJ? }, the 
two braiding operations n and T2 arc represented by 
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diag(l — i, 1 + i, 1 — i, 1 + i), 
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Based on these operations, we construct two composite 
braiding operations to demonstrate T1T2 7^ T2T1. The two 
composite braiding operations are chosen as S = TiTa, 
T = TiT\ with the properties 
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Suppose the two Majorana fermions arc initially pre- 
pared in the state (|0) + |1)) L ® (|0) + \V)) R /2, after 
the braiding operations ST and TS, the output states 
-i(|0) + |l)) £ 8) (|0) + |l)) fl /2 and -i(|0)-|l)) L ® 
(|0) + |l)) R /2 are orthogonal with each other. Thus we 
can demonstrate the non-Abelian nature of the Majo- 
rana fermions unambiguously by detecting the difference 
of the two orthogonal output states. 

In fact, this detection is readily available in experiment 
by transferring them into spin basis. In this basis, the 
two degenerate ground states of the topological lattice 
segment i are written as 
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ID 



V2 

(10) 

where |±) = (|t) ± The Dirac fermion 

and the fermion parity operators become c\ = (o~f — 

^m,nr=i^)/2 and P = nr=i CT l' where i = L,R, 
rrii is the total lattice sites in the topological seg- 
ment i. Using these expressions, one can demon- 
strate c+ 10> 4 = |l) <f c l \l) i = \0) i and P, !<)>,. fll),) = 
\0) i (— as in the JW fermions basis. By substitut- 
ing Eq. (10), the initial state of Majorana fermions is 
transformed into |_|_)® mL |_|_^® m «^ winch can be easily 
prepared by optical pumping, and the two output states 



become — i 1+) 



imp 



and i 



im L | + ^m B> which 

can be adiabatically rotated into -i ||)® mL an d 
• |^®m £ |^®m B by anothcr lascr driving. That is, we 
only need to distinguish the two orthogonal collective 
spin states \l)® mL and It)®" 11 *. For this purpose, one 
can apply one laser beam on the lattice segment L to 
couple the spin up state to an auxiliary excited state 
and detect the fluorescence of its emissions. The cases 
with and without fluorescence are corresponding to the 
state |t) and The distinct advantage here is that the 
state It)® 1 " 1, is a collective spin state which can provide 
a strong fluorescence to distinguish itself from the state 



|t)^ mi . Such property makes our proposal much more 
appealing for the detection the non-Abelian statistics of 
Majorana fermions. 

In summary, we have proposed an experimental scheme 
to create and manipulate Majorana fermions. In the 
scheme, the transverse field Ising model has been sim- 
ulated using cold atoms trapped in optical lattice. By 
tuning the intensity of the addressing lasers, Majorana 
fermions can be generated at the ends of the topological 
lattice segment. Such tuning can also allow them to be 
moved, fused and braided. Finally, we have shown how 
to detect the non-Abelian statistics of Majorana fermions 
by distinguishing two orthogonal output collective spin 
states immediately after performing two opposite order- 
ings of braiding operations. 
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